Inflation driven by causal heat flux 
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■ Abstract 

in , 

C*~> , We find a simple inflationary solution in an inhomogeneous spacetime with heat 

flux. The heat flux obeys a causal transport equation, and counteracts the infla- 
tionary decrease of energy density. At late times, the heat flux tends to zero and 

r**"* , the fluid approaches the equation of state p = —p. 
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1 Introduction 

Inflationary expansion arises when the effective pressure becomes sufficiently negative 
that its 'repulsive' contribution to gravity becomes dominant. This scenario occurs nat- 
urally in scalar-field models of the early universe. On the other hand, various papers 
have considered the possibility that inflation could occur in an 'ordinary' fluid if there is 
sufficient bulk viscous stress to drive the effective pressure negative (see, e.g., 0, |3| and 
references given there). Clearly, as pointed out in ||, the bulk viscous stress must exceed 
the local equilibrium pressure, so that such models are far from equilibrium. 

Bulk viscous stress is a scalar dissipative effect, and thus compatible with the spatial 
homogeneity and isotropy of Friedmann-Lemaitre-Roberston- Walker (FLRW) spacetime. 
The vector dissipative effect of heat flux is not compatible with FLRW symmetry. Here we 
present a 'toy' model to show that, in a simple inhomogeneous geometry, it is in principle 
possible for inflationary expansion to arise in a fluid with heat flux. Essentially, the heat 
flux acts to counter the decrease of energy density p due to expansion, while the pressure 
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p is steadily reduced (with asymptotic limit p — > — p). Unlike the viscous inflationary 
models, the fluid is not far from equilibrium at all spacetime events. 

For consistency, we use the causal transport equation of Israel and Stewart for heat 
flux, in preference to the often used Fourier equation (generalized to relativity), which is 
noncausal and has unstable equilibrium states ||. 



2 The basic equations 

Consider Modak's shear-free spherically symmetric model, given in comoving coordinates 
by! 

ds 2 = - [l + M(t)r 2 ] 2 dt 2 + a 2 (t) [dr 2 + r 2 (d6 2 + sin 2 Odcf)] , (1) 
where the fluid four- velocity is u a = [1 + Mr 2 ]~ 1 5 a Q. The fluid four-acceleration is 

2Mr 



pu>a 



1 + Mr 2 



5a 1 , 



and the volume expansion rate is 

= V a u 



3H 



1 + Mr 2 

where a dot denotes d/dt. The heat flux vector is 



H = 



^ 5 a i where q 2 (t,r) = q a q a 



(2) 



(3) 



(4) 



Note that q is a covariant scalar measure of the heat flux magnitude. The four-velocity 
u a is comoving with the particle frame, in which the total energy flux is the heat flux, 
since there is no particle flux relative to this frame. 
The Einstein field equations^ are satisfied if 



P 
P 

q 
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(1 + Mr 2 ) 2 
4MrO 



9 • - a 4M(1 + Mr 2 ) 

|Mr 6 - 2- + K — '- 

6 a a 2 



3(1 + Mr 2 )a 
Then the conservation equations 



(5) 
(6) 

(7) 



u a V a p +(p + p)Q + V a q a + a a q a = , 

(p + p)a a + hj {Vpp + tfV^p) + |6g Q = 



(9) 



are identically satisfied, where h a p = g a p + u a u$ projects orthogonal to u a , and g a p is the 
metric. In Modak's model, the heat flux is radially inward if the fluid is expanding, and 

1 We use units with 8irG = 1 = c 
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this is reflected in the fact that the energy density at each instant of time is a maximum 
at the centre of symmetry. 

The functions a(t) and M(t) are determined after appropriate thermodynamic equa- 
tions are imposed. We now require that the heat flux is governed by the causal transport 
equation of Israel and Stewart || 

Th aP u l V ry qp + q a = -n (h a/3 VpT + Ta a ) , (10) 



where T is the local equilibrium temperature, k (> 0) is the thermal conductivity and 
t is the relaxational time-scale which gives rise to the causal and stable behaviour of 
the theory. The noncausal Fourier transport equation has r = and reduces from an 
evolution equation to an algebraic constraint on the heat flux. Intuitively, one can see 
that in this case, the heat flux is instantaneously brought to zero when the temperature 
gradient and acceleration are 'switched off'. 

For the line element (p, equation (jT^) becomes 



rq+ (l + Mr 2 )g = -- [(l + Mr 2 )^' , (11) 
where a prime denotes 8/ dr. 



3 A simple inflationary model 



To close the system of equations, one needs further thermodynamic information about r, 
k, p, p and T. In ||, this is done for a radiation-dominated model where heat transport 
arises from radiative transfer. Here we follow a different approach, since we are interested 
in demonstrating the theoretical possibility of inflation driven by causal heat flux. First 
we suppose that M = M is a positive constant. Then on each comoving sphere r = 
constant of fluid particles, inflationary expansion is characterized by a > 0. A particular 
case of inflationary expansion is H — Hq where Hq is a positive constant, and then 
a = aoexp(H t). Substituting into equations (0)-(0), we find 
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V 
Q 



3H* 



(1 



M r 2 ) 2 ' 
4M 



a 2 (1 + M r 2 ; 
4M H r 
a (1 + M r 2 ) 



e- 2Hot -p, 



-H t 



(12) 
(13) 
(14) 



Note that the magnitude of the heat flux is a maximum at r = 1/ ^3M , falling to zero 



at the centre and as r — >• oo. The causal transport equation (|TT| ) becomes 



AM H r [(l + M r 2 ) - H r] = k (l + M r 2 ) 2 [(l + M r 2 ) T]' . (15) 

We have yet to specify k, t and equations of state involving the temperature. However, 
we will avoid the difficult issue of trying to introduce a microscopically motivated model, 
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since our primary aim is only to show the possibility of consistent solutions. In this spirit, 
we will satisfy the causal transport condition ( |i"5| ) by taking 



T 



(l + M r 2 ) if" 1 , (16) 

T = U ® (17) 
1 + M r 2 ' { ' 

where U is an arbitrary positive function. Equation ([H]) implies that the relaxation time 
is the same as the local expansion time, which is not unreasonable. The radially inward 
heat flux counters the cooling brought about by expansion. However, one might expect 
that the inflationary cooling will not be balanced or overcome by dissipative heating, and 
then one would choose U < 0. 

For any choice of U, it follows from ([L7|) that T decreases radially outward, in apparent 
contradiction to the radially inward direction of the heat flux. However, in relativity, the 
inertia of heat energy gives rise to an accelerative contribution to the heat flux. In our 
model, the latter dominates the temperature gradient. Note that, whatever the (positive) 
choice of U and k, the second law of thermodynamics will be satisfied, since it is built 
into the causal theory p. 

The model given by equations (|P2|)-(|T4|), ( |T6| ) and ( |ITD describes inhomogeneous infla- 
tion driven by causal heat flux. The inflationary expansion rapidly reduces the pressure, 
as shown by equation (|13|) , but the heat flux counters the tendency for the energy den- 
sity to decrease, and equation (|12|) shows that the energy density is a comoving constant 
{u a V a p = 0). At late times, the heat flux tends to zero and the fluid approaches a 
de Sitter-like equilibrium with p = —p. The deviation of the fluid from equilibrium is 
measured by the covariant dimensionless ratio 

M = (18 ) 

p \Sa H / 

At late times, this ratio rapidly becomes small. However for any fixed finite time, the 
ratio grows with radius, reflecting the fact that \q\ decays less rapidly than p. Thus the 
fluid is only close to equilibrium near the centre and at late times. By contrast, in bulk 
viscous inflation the fluid is far from equilibrium at all times and positions. 

We have shown that in principle it is possible for causal heat flux to drive inflation. 
This approach could be the basis for a model of inflationary expansion of bubbles in the 
early universe, driven by causal energy transport processes. Our model is of course highly 
simplified, and lacks some of the physical properties that would be expected of a more 
realistic model. For example, there is no mechanism in our simple model for achieving an 
exit from inflation, and the model inflates for all time. Furthermore, the heat flux does 
not homogenize the universe, as one might intuitively expect, so that the metric does 
not tend to the FLRW metric asymptotically. This is a direct consequence of the simple 
choice of M(t) as a constant. However, in principle it should be possible to overcome 
these drawbacks with a more sophisticated and complicated model. 
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